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Abstract 

The comparison of B-meson production cross sections as the results of PYTHIA 
code and Quark Gluon String Model (QGSM) is carried out for energies of proton 
colliders: Sp(p)S, Tevatron and LHC. Model predictions are based on the idea of 
supercritical Pomeron exchanges with the phenomenological intercept Ap(0)=0,3 
for heavy quark production. Transverse momentum spectra of B-mesons are also 
compared. It is shown that the cross sections calculated with PYTHIA using CTEQ 
structure functions are in a contradiction with the asymptotical estimation of BB 
production cross sections in QGSM. Asymmetries between the spectra of B° and 
BO mesons are also contradicting. The reasons of the difference are discussed. 
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1 Introduction 



We can't say that the complete knowledge on beauty quark pair production 
is obtained now, since the data we have on bb production cross sections are 
not yet sufficient. In order to monitor the model ideas on the phenomenon 
it seems useful to revise once in a while the collected data. Recently the 
results of several experiments [jl], g, |3j carried out at two energies of colliding 
protons, 630 GeV and 1.8 TeV, are represented in the literature. 

In this article the two of the models are compared: on the one hand 
- the phenomenological Quark-Gluon String Model ||], based on idea of 
hadronic amplitude duality and on the theory of supercritical Pomeron , 
on the other hand - the wide-spread Monte Carlo code PYTHIA j5j] , which 
includes the results of QCD perturbative diagram calculations. 

The production cross sections arising with energy is a fact which was 
widely discussed in recent studies [§, PPfl . The theory of supercritical 
Pomeron postulates the rising as s Ap (°\ where Ap(0) = op(0) — 1, ap(0) 
is the intercept of Regge trajectory of Pomeron. 

The energy behavior of the production cross section in QCD perturba- 
tive approach are provided by the choise of gluon structure functions of 
interacting hadrons. Most of those functions, which are accepted now for 
MC simulations of high energy collisions, are built to approximate the re- 
cent data from HERA measured up to x = 1CT 4 . It should be noticed, that 
all known gluon structure functions which satisfied this recent data, can be 
taken for the modelling of bb production at LHC, because the value 10~ 4 
is the very region of x attributed for B meson production at 14 TeV due to 
the following estimation: 2rap/(14Tey) ~ 10~ 4 . One of those appropri- 
ate gluon distributions is CTEQ structure function, which is involved into 
PYTHIA code as default distribution. 

2 Parameters defining the B-meson production cross 
sections in QGSM. 

The major parameter of QGSM which defines the cross section dependence 
on energy is Ap(0) e //, which have to be less than the same value for one 
Pomeron exchange, because multipomeron diagrams or branches should 
be taken into account in the calculations. This parameter depends on the 
mean value of transverse momenta transmitted in the process. Therefore 
the energy dependences of different mass particle production cross section 
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must be led by different Ap(0) e //. 
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Figure 1: The A(Q 2 ) dependence obtained in HI |7[] and one Pomeron exchange approx- 
imation (the solid line). 

The Ap(0)(Q 2 )-dependence can be illustrated well with the data ob- 
tained in HI experiment at HERA jjfl. The Pomeron exchange plays an 
important role in electron-proton collisions too. The multyhadron pro- 
duction takes place in this process due to the Pomeron exchange between 
photon and beam proton. The A(Q 2 ) data are shown in Fig. , where each 
point was extracted from measured F2(x,Q 2 ) function by the approxima- 
tion with the simple dependence: ~ x^ A ^ K 

It should be noticed that the theoretical curve was calculated [|7j un- 
der the assumption of one Pomeron exchange which doesn't include the 
branches as against the proton-proton interaction. The QGSM scheme for 
heavy meson production have to be similar to the one Pomeron exchange 
pattern due to the energy conservation reasons. So we can take for the A e ff 
the value 0,3 corresponding to the HI data approximation at Q 2 = (2m_e) 2 . 
It's worth beeng stressed here, that this value differs from A e //=0.12 which 
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was chosen for light quark meson production in early papers |TTj . 

Another Regge trajectory parameter is important for model description 
of inclusive cross sections of B meson pair production, it's ax(0), the inter- 
cept of 66-trajectory. It provides the increasing cross section at the energy 
region close after the B pair production threshould. There are various opin- 
ions about the value of this parameter. From the point of view of QGSM 
it may vary in the range of -16-1-0 |n|. The other athors prefere to take 



a T (0) = -9 

The value ax(0) = — 16 will be taken here for to estimate the upper limit 
of growing cross section when it increases rapidly after the threshould. The 
parameter discussed above exists in the functions of fragmentation of qurk- 
gluon strings into each sort of B-mesons. Those functions are written in 
QGSM according to the rules fulfilled by the Regge asympthotics ||. 

For example, the function for d quark string fragmentation into B + 
contains the following factors: 

Vf(z) = ^(1 - s)-M0)«r(0)+A (1 + a B z 2^ 

where is the density parameter for fragmentation of quark-gluon 
string into B-mesons. The af is the parameter of string fragmentation 
asymmetry introduced in [jTUj] to provide a transition between probabilities 



of the B production at z^O and z— »l.The value af can be of the order 10 
and actually doesn't impact on the value of B production cross section at 
energies higher than 1.8 TeV. 

The calculations of p± spectra of produced hadrons are also available in 
the framework of QGSM, as it had been done already in for 7r-mesons. 



The spectra can be described up to the momenta of order few GeV/c 
in this substantially nonperturbative model . The distributions were of 
exponential character at low p± in this approach. Therefore the transverse 
momenta distributions for heavy flavor particles was not elaborated in this 
model. 



3 PYTHIA machinery 

The version PYTHIA 5.7 was taken to calculate the spectra of B-mesons 
at three energies of colliding protons: 630 GeV, 1,8 TeV and 14 TeV. The 
CTEQ gluon structure function [Jnj are used in this version to discribe the 



increasing cross sections. On the one hand the process of production of 
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such heavy quarks as b is good enough for beeing described by perturbative 
QCD diagram with gluon-gluon fusion. On the other hand, more and more 
low x gluons are involved into this process at energy rising and cross section 
becomes dependent on the accuracy of gluon structure function measured 
at low x. 

As it was mentioned in Introduction, we have precise data on due 
to HERA experiments up to x ~ 1CT 4 , what is enough for the calculation 
of bb production at LHC energies. Such a way CTEQ structure functions 
have to provide the right description of increasing cross sections of bb pair 
production. 

However, b quarks can be obtained not only in gluon fusion process. 
Two additional ways exist to produce bb pair, they are: gluon splitting 
99 ~~ * 99 1 where gluons gives bb pair in the next- to-leading order of correc- 
tions, and heavy flavor exitation Qig — > Qig. In PYTHIA this subprocesses 
are taken into account with massless matrix elements. It is a problem how 
to sum the resulting distributions from such different deposits. 
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Figure 2: Transverse momenta distributions of B-mesons fitted with PYTHIA. 
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It makes the ^-spectra at 1.8 Tev comparable with the data obtained 
in CDF experiment ( see Fig.). At the same time there is not good de- 
scription of UA1 data . It looks like the p± spectra were increased with 
additional fractions only by a factor and there is not any difference between 
the patterns of spectra for different subprocesses. It leads to rather flat 
form of transverse momentum distributions in PYTHIA and to small total 
cross section of B production at various energies. 

4 Comparison of cross section energy dependences 

The resulting energy dependences of production cross section are shown 
on Fig. for PYTHIA program and for QGS model as well. 
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Figure 3: Energy dependences of B-meson cross section. 

As it was mentioned above the QGSM curve is highest as it is possible 
in this model after the normalisation to the CDF cross section. But latter 
was obtained due to PYTHIA (see Fig. ).Thus the point where both curves 
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are crossing is rather conventional and depends on the form of ^-spectra 
at low transverse momenta accepted in PYTHIA. 



5 PYTHIA and QGSM predictions for the asymme- 
try between B°/B° spectra 

It would be interesting to consider the leading effect in the spectra of B- 
mesons at various energies. The valuable asymmetry between B°- and B°- 
meson spectra at LHC energy will be important for CP violation measure- 
ments. The recent prediction of the y dependence of such leading/nonleading 
asymmetry |T3j provided with PYTHIA simulations gives zero value of A(y) 
in wide range of y at the central region (see Fig. ). 



0.5 




Figure 4: B°/B° asymmetries at LHC energy given by PYTHIA [|T^] and by QGS model; 
the value A e //=0.3 correspods to the production of the mass of B-meson. 

The A(y) dependence in fragmentation region xp — ► 1 contradicts with 
all similar asymmetry measurements for D-meson spectra [[15], [16], [T7|, [T3j . 
The intersection of inclusive spectra of different type of B-mesons gives 
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the asymmetry passing trough zero at some yo, while the measured spectra 
of leading particles are always higher then nonleading particle spectra, so 
the asymmetry is positive value. In opposite to PYTHIA predictions, the 
asymmetry calculated in the framework of QGSM is rising function up to 
xp — > 1. This behavior is usually peculiar for the string approach because 
of so called "beam drag" effect. The valuable asymmetry in central region 
given in QGSM prediction jTBj is not enough small for not to be taken into 
account at CP violation measurements. It looks important to consider 
both these predictions in details and to discuss the probability of nonzero 
asymmetry in the production spectra at LHC energy. 



6 Conclusions 



We have compared two approaches for the understanding of the heavy fla- 
vored particle production: one of them is mostly perturbative and another 
one is nonperturbative at all. This comparison shows that some different 
suggestion has to be done for low transverse momenta distributions of B- 
mesons to put into agreement the both model predictions at LHC energy. 
The contradicting dependences for B / B° asymmetry in the B meson pro- 
duction spectra might be important for the CP violation measurements. 
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